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Abstract 

A search for an excess of muon-neutrinos from neutralino annihilations in the Sun 
has been performed with the AMANDA-II neutrino detector using data collected in 
143.7 days of live-time in 2001. No excess over the expected atmospheric neutrino 
background has been observed. An upper limit at 90% confidence level has been 
obtained on the annihilation rate of captured neutralinos in the Sun, as well as the 
corresponding muon flux limit at the Earth, both as functions of the neutralino 
mass in the range 100 GeV-5000 GeV. 

Key words: Dark matter. Neutrino telescopes, AMANDA, neutralino. 
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PACS: 95.35.+d, 95.30.Cq, ll.30.Pb 



1 Introduction 



The Minimal Supersymmetric extension to the Standard Model of particle 
physics (MSSM) [1] provides a promising dark matter candidate in the hghtest 
neutralino, x°, a linear combination of the supersymmetric partners of the 
electroweak neutral gauge- and Higgs bosons. Assuming R-parity conservation, 
the neutralino is stable. From accelerator searches and relic density constraints 
from WMAP data, a lower limit on the mass of the MSSM neutralino can be 
derived [2]. Typical lower limits for m~o from such studies are about 20 GeV, 
depending on the values chosen for tan/3. In models where the pseudo-scalar 
Higgs boson A is assumed to be light, < 200 GeV, the neutralino mass 
can be as low as 6 GeV. Theoretical arguments based on the requirement 
of unitarity set an upper limit on m~o of 340 TeV [3]. Within these limits, 
the allowed parameter space of minimal super-symmetry can be exploited to 
build realistic models which provide relic neutralino densities of cosmological 
interest to address the dark matter problem. 

Relic neutralinos in the galactic halo can become gravitationally bound in 
orbits in the solar system by losing energy through elastic scattering with 
matter. They may finally be trapped and accumulate inside celestial bodies 
like the Sun [4,5,6], annihilate, and produce a neutrino flux from the decays 
of the annihilation products [7] . 

We have previously published a search for neutralino dark matter accumu- 
lated in the Earth using the 10-string detector AMANDA-BIO and data from 
1997 [8]. In this letter we present the first search for neutralino dark matter 
accumulated in the Sun with the extended AMANDA-II detector. Due to the 
position of the Sun at the South Pole, reaching no more than 23.5° below the 
horizon, the separation of a potential neutralino-induced neutrino flux from 
the Sun from the atmospheric muon flux is a challenge that we could only 
address because of the increased sensitivity of the AMANDA-II detector to 
horizontal tracks. 



Email address: cpli@tsl.uu.se (C. P. de los Heros). 
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2 The AMANDA-II detector 



AMANDA-II consists of an array of 677 optical modules deployed on 19 verti- 
cal strings at depths between 1500 m and 2000 m in the South Pole ice cap. The 
strings are arranged in three approximately concentric circles of 60 m, 120 m 
and 200 m diameter respectively. An optical module consists of a photomul- 
tiplier tube housed in a glass pressure vessel. Muons from charged-current 
neutrino interactions near the array are detected by the Cherenkov light they 
produce when traversing the ice. During 2001 the detector was triggered when 
at least 24 modules were hit within a time window of 2.5 yus. The trigger rate 
was 70 Hz. The relative timing of the photons reaching the optical modules 
allows the reconstruction of the muon track. A more detailed description of 
the reconstruction techniques used in AMANDA is given in Ref. [9] . 



3 Signal and background simulations 

The simulation of the neutralino-induced neutrino signal was performed using 
the DARKSUSY program [10] for a sample of neutralino masses (100, 250, 500, 
1000, 3000 and 5000 GcV). Neutralinos can annihilate pair-wise to, e.g., 
qq, W+W-, Z°ZO, Rl^^l Z^B.^^ and W±H^, and neutrinos are produced 
in the decays of the annihilation products. Neutrinos produced in quark jets 
(from e.g. bb) or from Higgs bosons typically have lower energy than those 
produced from decays of r leptons or gauge bosons. Two annihilation channels 
were considered in this paper for each mass tested, XiXi ~^ bb and XiXi ~^ 
W"'"W~, which are referred to as soft and hard, respectively, in the rest of 
the paper. This choice covers the range of neutrino energies which would be 
detectable with AMANDA for typical MSSM models. The simulated angular 
range was restricted to zenith angles between 90° (horizontal) and 113°, with 
the generated number of events as a function of angle weighted by the time 
the Sun spends at each declination. The neutrino-nucleon interactions in the 
ice around the detector producing the detectable muon flux were simulated 
with PYTHIA [11] using the CTEQ3 [12] parametrization of the nucleon structure 
functions. 

The background for this search arises from up-going atmospheric neutrinos 
and mis-reconstructed downward-going atmospheric muons. We have simu- 
lated the atmospheric neutrino flux according to [13]. The simulation includes 
neutrino propagation through the Earth, taking into account the Earth den- 
sity profile, neutrino absorption and neutral current scattering. A sample of 
1.9x10"^ atmospheric neutrinos with energies between 10 GeV and 10 PeV and 
zenith angles between 70° and 180° (vertically up-going) has been simulated. 
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Fig. 1. Relative cut efficiencies for one of the neutralino masses and annihilation 
channels considered in this paper, 500 GeV hard. The efficiency is set to 1 at trigger 
level, except for the atmospheric neutrino curve, that has been scaled relative to 
the number of events expected in the live-time of the analysis. The axis scale on the 
right gives the absolute normalization for the data and the atmospheric neutrino 
curves. 



The simulation of atmospheric muons was based on the CORSIKA air shower 
generator [14] using the South Pole atmosphere parameters and protons as 
primaries. We have simulated 1.6 x 10^ interactions, distributed isotropi- 
cally with zenith angles between 0° and 85°, and with proton energies, Ep, 
between 600 GeV and 10^^ GeV, assuming a differential energy distribution 
oc Ep^-"^ [15]. The sample corresponds to about 32.5 days of equivalent detector 
live-time. 



Muons were propagated from the production point to the detector taking 
into account energy losses by bremsstrahlung, pair production, photo-nuclear 
interactions and 5-ray production as implemented in the code MMC [16]. 
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4 Data analysis 



The data set used comprises 8.7x10^ events collected in 143.7 days of effective 
live-time between the end of March and the end of September 2001. The data 
and Monte Carlo events have been reconstructed with a fast first-guess track 
finding algorithm which is used as a seed for an iterative maximum-likelihood 
reconstruction. The reconstructed events were processed through a series of 
filters in order to reduce the atmospheric muon background and retain as 
much of a potential signal as possible. The variables used in the analysis are 
the same for the different neutralino masses, but the cuts have been optimized 
independently for each mass in order to exploit the differences in the resulting 
muon energy spectra. The event selection was optimized in a 'blind' way, by 
exphcitly excluding the time tag of the events (and thus the direction of the 
Sun) in the development of the selection criteria (see [17] for details of this 
analysis) . 

Firstly, events were selected according to the quality of the track reconstruc- 
tion and their zenith angle. The quality criteria were based on 'direct hits', i.e. 
hits with time residuals, t^es-i in the interval [—15ns, 25ns] relative to the ex- 
pected arrival time of an unscattered photon from the fitted track. Cuts were 
made on the number of direct hits, the maximum distance along the track 
direction between two such hits and the number of strings with direct hits. 
This selection (cut levels 1 to 3 in figure 1) resulted in a rejection of 99.97% 
of the data (mostly downgoing muons) , while the simulated background from 
atmospheric neutrinos was reduced to 30%. The simulated signal was reduced 
between 46% and 54% (depending on the neutralino mass). 

The next step consisted on a cut on a neural network (NN) output. For each 
neutralino mass, a NN was optimized for the simulated signal and a back- 
ground sample consisting of simulated atmospheric muons. No information 
about the characteristics of atmospheric neutrinos was included in the training 
of the network. Eight variables were used as an input to the NN to characterize 
the events: the number of hit channels, the difference in reconstructed angle 
between the fast first-guess fit and the maximum likelihood fit, the number 
of 'late' {tres >150 ns) hits, the track length of the first-guess fit, the length 
spanned by all the hits projected onto the direction of the final reconstruc- 
tion, the likelihood of the fit, the center of gravity of the hits and the vertical 
distance between the deepest and shallowest hit modules. The MLPfit pack- 
age was used [18]. The network performance was optimized by varying the 
network architecture, the size of the training samples and minimization algo- 
rithms. An architecture of 8-40-1 with a hybrid linear minimization was found 
to give maximal efficiency. 

The final event selection consisted of the following series of cuts: the number 
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Table 1 

For each ncutralino mass and annihilation channel the table shows: The number of 
data events, N^ata, the number of expected atmospheric neutrinos in the declination 
band of the Sun, N^^^,, the search bin, the final number of data events and the 
estimated background in the search bin, N*^^ and N*^ respectively and the muon 
effective volume, Vca. The last three columns show the 90% CL upper limit on the 
expected signal, figo, and the corresponding 90% CL limits on the annihilation rate 
at the center of the Sun and on the muon flux at the Earth, Fa and The limits 
include systematic uncertainties. 'Soft' and 'hard' refer to annihilation into bb and 
W^W~ respectively. 
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of hit channels, the number of 'early' hits (t^^s < -15 ns), the number of hits 
used in the first-guess fit and the distance between the first and last direct hits 
projected on the track direction. For this last cut, the direct hit definition used 
was -15 ns < t,^, <75ns. Figure 1 shows an example of the relative efficiency 
for signal, background and data as a function of cut level, for one of the 
neutralino masses. Cut level 4 corresponds to the cut on the NN output, and 
the last level to the addition of the cuts just described, shown individually 
as the intermediate points between level 4 and 5 in the plot. Figure 2 shows 
the distribution of the cosine of the space angle, fl, between the reconstructed 
tracks and the direction of the Sun (after unblinding). The distribution is in 
agreement, both in shape and normalization, with the expectation from the 
atmospheric neutrino background, within the 25% overall systematics of the 
atmospheric neutrino flux [19] (indicated by the shaded band). 
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Fig. 2. Cosine of the space angle between the reconstructed track and the Sun 
position at final cut level for data and atmospheric neutrinos. The case shown 
is for the cut optimization for 500 GeV neutralinos (annihilation into 
The shaded area shows the one standard deviation systematic uncertainty in the 
expected atmospheric neutrino distribution. 

5 Systematic uncertainties 



There are several systematic uncertainties in the detector effective volume 
that should be taken into account when calculating limits. The uncertainty 
in the optical module sensitivity contributes about 14%, while the detector 
calibration and the hardware simulation contribute about 5%. The largest 
contribution comes from the imprecise knowledge of the layered structure and 
the impurity contents of the ice, and is less than about 30% (with the exception 
of the softest spectrum, annihilation of 100 GeV neutralinos into bb where this 
uncertainty reaches about 50%). Uncertainties on the neutrino- nucleon cross- 
sections and in the simulation of muon energy loss contribute of the order of 
5%. The overall systematic uncertainties in the detector effective volume lie 
between 20% and about 60% depending on the neutralino mass model being 
tested (see [17] for details). The uncertainties have been incorporated in the 
limits using the method described in [20] with the unified Feldman- Cousins 
ordering scheme [21]. 
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6 Results 



Different angular bins around tfie position of tlie Sun for each neutralino mass 
and annihilation channel were used in the search for a signal. The bins were 
chosen as containing 90% of the signal, and they range from 5.5° for the harder 
neutrino spectra to 26° for the softest spectrum (100 GeV neutralinos anni- 
hilating into bb) . The off-source data in the dechnation band of the Sun were 
used to estimate the expected background at the Sun position. This procedure 
eliminates any effects of uncertainties in the background simulation and in the 
total normalization of the atmospheric neutrino flux. Table 1 summarizes the 
relevant numbers. 

No evidence of a statistically significant excess of events from the direction of 
the Sun was found. Therefore, the number of observed events from that di- 
rection, together with the Poisson mean of the background expectation, were 
used to set limits on the neutrino flux at the Earth from neutralino annihi- 
lations in the Sun. We have followed the same prescription as in [8], starting 
from the directly measurable quantity, the neutrino-to-muon conversion rate, 
Ti,^n = iigo/ (V^ff et), where /xgo is the 90% CL signal upper limit, t is the live- 
time, e the reconstruction efficiency and V^g the effective volume of the detector 
at final cut level. The limit on the annihilation rate of neutralinos in the Sun, 
Fa, is proportional to T^^^, where the proportionality coefficient takes into 
account the production of muons through the neutrino- nucleon cross-section 
weighted by the different branching ratios of the XiXi annihilation process 
and the corresponding neutrino energy spectra. Prom Pa we reach the limit 
on the muon flux at the Earth by taking into account that 



where Dq is the distance to the Sun and dN^/dE^ is the muon flux produced 
from the neutralino annihilations, including all MSSM model dependencies 
and interaction kinematics, ii^thr is here an arbitrary threshold, that can be 
used to convert the measured flux limit (which is obtained for the actual de- 
tector energy threshold) to any other threshold. The last three columns of 
table 1 show the Feldman-Cousins 90% upper limit on the expected signal 
and the corresponding 90% CL limits on the annihilation rate at the center 
of the Sun and on the muon flux at the Earth. All numbers include system- 
atic uncertainties. 
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Neutralino Mass (GeV) 



Fig. 3. Upper limits on the muon flux from the Sun, from neutrahno annihi- 
lations into W~^W~ (hard channel) as a function of neutralino mass. The muon 
energy threshold has been extrapolated to a common value of 1 GeV. The dots 
show the models disfavored by recent direct search results from CDMS. 

7 Discussion 



Figure 3 shows the AMANDA limits on the muon flux at the Earth from 
neutralino annihilations into W^W~ (hard channel) in the Sun as a function 
of neutralino mass, compared to the results of Baksan [22], MACRO [23] 
and Super-K [24] and theoretical predictions based on DARKSUSY. The limits 
have been rescaled to a common muon threshold of 1 GeV using the known 
energy spectrum of the neutralinos. The sparseness of the arrangement of the 
AMANDA strings (not less than 25 m between nearest neighbours) makes 
the detector less effective at lower neutrino energies, and leads to a worsening 
of the limit for neutralino masses below about 500 GeV. Above this mass. 
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AMANDA, with the 144 days of live-time used for the present analysis, is 
already competitive with experiments with much longer exposure times (3.6 
years of collected data for MACRO and 4.6 years for Super-K for example). 

The symbols in the figure represent the predictions of DARKSUSY for given 
combinations of SUSY parameters. The models shown are only those which 
give a relic density within the current limits of allowed dark matter density 
in the universe, 0.05< Qh!^ < 0.2. The dots represent parameter combinations 
that are disfavored by the latest results from direct experiments setting lim- 
its on the neutralino-nucleon cross-section [25] . Although comparison between 
direct and indirect searches is not straightforward, the figure shows that cur- 
rent results from indirect searches are competitive in the high neutralino mass 
region. 

Neutrino oscillations in the atmosphere (not included in the simulations) 
would tend to reduce the muon-neutrino flux and affect the total normal- 
ization of the background. Since we have used the measured neutrino flux 
off-source as an estimation of the background at the source position, the effect 
of oscillations is inherently included in the analysis. The atmospheric neutrino 
simulations were used as a confirmation that the analysis cuts selected at- 
mospheric neutrinos and rejected the down-going atmospheric muons to the 
desired level before the unbhnding of the data set. 

The DARKSUSY simulations did not include oscillations of the neutrinos on 
their way to the Earth. The effect of oscillations of neutrinos produced in the 
Sun on the final fiux observed at the Earth has been discussed in [26] (and 
in [27] in a slightly different context). Detailed calculations in [26] found that 
the resulting muon flux at the Earth from low mass {m~o < rrit) MSSM neu- 
tralino annihilations can be considerably affected by Ut- oscillations in 
some scenarios. This is the case for MSSM models with important annihilation 
branching ratios into t~^t~ , typically Brr > 0.1. In such cases the muon flux 
at the Earth can be increased by a factor of about 2 to 4 with respect to the 
no-oscillation case. For higher neutralino masses, where the main annihilation 
channel can be through gauge bosons or heavy quarks, the effect of Vt i/^ 
oscillations is negligible, oscillations into z/^ being compensated by i/^ oscillat- 
ing into z/^. Flavor mixing with z/g both from z/^ and z/^ is suppressed by the 
small value of ISm^^ s-nd the mixing angle sin^ 2^^3i respectively, and does not 
play a role when considering the flux from the Sun. Given the considera- 
tions above, the flux and annihilation rate limits presented in this letter for 
the case of annihilations into W~^W~ would be practically unaffected if one 
were to consider oscillations. The results for the bb annihilation channel would 
also be unaffected for high neutralino masses {rn~o > 200 GeV) and would 
worsen by a few percent for neutralino masses below 200 GeV. 
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